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Abstract — The constraints on total neutrino mass and effective number of neutrino species based on 
CMB anisotropy power spectrum, Hubble constant, baryon acoustic oscillations and galaxy cluster mass 
function data are presented. It is shown that the discrepancies between various cosmological data in 
Hubble constant and density fluctuation amplitude, measured in standard ACDM cosmological model, can 
be eliminated if more than standard effective number of neutrino species and non-zero total neutrino mass 
are considered. This extension of ACDM model appears to be ss 3<r significant when all cosmological data 
are used. The model with approximately one additional neutrino type, N e s ~ 4, and with non-zero total 
neutrino mass, Em„ ~ 0.5 eV, provide the best fit to the data. In the model with only one massive neutrino 
the upper limits on neutrino mass are slightly relaxed. It is shown that these deviations from ACDM model 
appear mainly due to the usage of recent data on the observations of baryon acoustic oscillations. Larger 
than standard number of neutrino species is measured mainly due to the comparison of the BAO data with 
direct measurements of Hubble constant, which was already noticed earlier. As it is shown below, the data 
on galaxy cluster mass function in this case give the measurement of non-zero neutrino mass. 
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INTRODUCTION 

The existence of neutrino with masses in the eV 
mass range would produce a suppression of density 
fluctuations on scales below the horizon when these 
neutr inos become non-rclativistic (e.g., iHu et all 
Il998 ). The amount of this suppression can be mea- 



sured from the comparison of density fluctuations am- 
plitude measured at early epoch from CMB anisotropy 
power spectrum and at recent epoch, from, e.g., galaxy 
cluster mass function data, which finally could give the 
measurement of total neutrino mass. 

From cosmological observations the number of neu- 
trino species could also be measured. Any type of neu- 
trino, which was thermalized in early Universe, make 
its own contribution into the energy density of rel- 
ativistic matter before equipartition. The change of 
this energy density result in a change in Universe ex- 
pansion rate and therefore in changes of the size of 
sound horizon and photon diffusion scale (Silk damp- 
ing scale). 

All these quantities can be measured using the data 
from various cosmological observations. For example, 
the information on density fluctuation amplitude at 
early epoch, size of sound horizon, Silk damping scale 
is contained in CMB anisotropy power spectrum which 



is now accurately m easured in various experiments 
such as WMAP 



■g- : 



Larson et al 



Telescope (jKeisler et all 1201 lUStorv et al 



20111). South Pole 
20121) At- 

acama cosmology telescope (jDunklev et all 120111 ). in 
very near future the results of CMB measurements ob- 
tained wit h Planck space observatory should also be 
published (jPlanck Collaboration , 120111) . 

The linear density fluctuations amplitude at recent 
epoch, which is usually described using a% parame- 
ter, may be measured by means of various methods. 
One of the most accurate measurement of this quan- 
tity comes fro m the data on galaxy clu st er mass func- 
tion (see, e.g. , Vikhlinin et all l2009al fbl; iMantz et al.l 



2010at iBenson et all l201ll: iBurenin. Vikhlininl . 120121) . 



The information on the size of sound horizon is 
contained also in the baryon acoust ic osci llations 



data (BAO, see, e.g., Percival et al. , 2010l) . In 



addition, t he direct m e asurem ents of Hubble con- 
stant (e.g., iRiess et ah . 2011 ) should also be used 
in order to eliminate degeneracies between var- 
ious cosmological parameters. The combination 
of all the data discussed above allowed to ob- 
tain constraints on total neutrino mass and mini- 



ber of neutrino species (e.g., Vi 


dilinin et al.. 2009b 


Mantz et al. 2010bUKeisler et al. 


,2011 


: Benson et al. 


201lt Burenin. Vikhlininl. l2012h. 
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In addition to the data used in these works, substan- 



1 



R. A. Burcnin 



tially improved new data on BAO ob s ervations were 



2011 



publi s hed recently (Blake e t al.l. 120111: iBeutler et al 



IPadmanabhan et al.l . 2012t Anderson et al 



20121 ). It was noticed that the distance scale measured 
using these data in assumption of standard ACDM 
model turns out to be in some tension with the re- 
sults of direct Hubble constant measurements and that 
this discrepancy can be eliminated with the assump- 
tion of larger than standard effective neut rino species 



(jMehta et all . l2012t lAnderson et all 120121) . With new 



additional data on Hubble constant measurements 



taken in account (jFreedman et all 12012!) , this discrep 



ancy have increased. 

As it is shown below, when these new BAO data are 
used, there is also a discrepancy in og measurements 
inferred from the data on CMB anisotropy power spec- 
trum in assumption of standard ACDM model and 
obtai ned using the g a laxy cl uster mass function data 
from IVikhlinin et~al1 (|2009allbh . This discrepancy can 
be interpreted as a suppression of density fluctuations 
due to the existence of neutrinos with non-zero total 
mass. It is shown that the best agreement with all cos- 
mological data discussed above is achieved in model 
with non-zero neutrino mass and additional neutrino 
species. 

The cosmological parameters constraints were cal- 
culated from the simulations of Monte-Carlo Markov 
Chains, which were done using CosmoMC software 



(jLewis. Bridld . 120021) . version of Jan. 2012. In all Fig- 



ures below the contours at 68% and 95% confidence 
levels are shown. All numerical values of confidence 
intervals are given at 68% confidence level. 

COSMOLOGICAL DATA 

Galaxy clusters 

For our work the data on galaxy cluster mass func- 
tion measurements were ta ken without any changes 
from Vikhlinin et al 1 (|2009al lbh. In this work a sample 
of 86 massive galaxy clusters with masses measured 
with about 10% accurac y using Chandra observations 
( Vikhlinin et al ], l2009ah was used. Distant clusters, 
located atzw 0.4-0.9, were selected from 400d X-ray 
galaxy cluster survey , based on ROSAT pointing data 



(Bur cnin et all 12007 ). Clusters in local Universe were 



selected using ROSAT all sky survey (see details in 



Vikhlinin et all 12009a). 



Likelihood functions for this cosmological datasct 
are available at WW W 1 (see also details in 
Burenin. Vikhlininl 2012 ). Systematic uncertainties 
are not included in these likelihood functions. These 



uncerta inties are discussed in detail in IVikhlinin et al.l 
(|2009bl) and ca n be t a ken in account se p aratel y 
(jvikhlinin et~aH l2009bl: iBurenin. Vikhlininl . |2012| ). 



which is done below for all measurements where these 
data are used. This dataset is designated below as CL. 

CMB power spectrum 

For our work we used the d ata of 7-year obser 



Komatsu et al 



2011 



vations of WMAP o bservatory (jLarson et al 

' 20111) . Likelihoods were calculated 



using the software taken from archive 2 , version J^.l. 
Also we used the data at small angular scale CMB 
anisotropy obtained with South Pole Telescope (SPT, 
Keisler et al.l . [2oTTh . When using these data, the con- 
tribution of Sunyaev-Zcldovich effect and also the con- 
tributions of poisson and clustered sources were taken 
into account accord ing to the prescription of §4.1 in 
Keisler et al.l (|2011l) . These data taken together are 



designated below as CMB. 

Hubble constant 

In our work we used the Hubble constant mea- 
surement obtained using improved calibration of su- 
pernovae type la ab solute magn i tudes, H = 73.8 ± 



2.4 km s^Mpc -1 (|Riess et all . l201ll ). The error 



here includes both statistical and systematic uncer- 
tainties. In addition to this measurement, the re- 
sults of Carnegie Hubble Project were published re- 
cently, where the new Cepheid distance scale cal- 
ibration, obtained using the data of Spitzer Space 
Telescope, was applied to the Hubble Space Tele- 
scope Key Project data and the measurement Ho = 
74.3±1.5(stat.)±2.1(sy s.) km s^Mpc" 1 was obtained 
(iFreedman et all 120121) 



As compared to iRiess et al.l (|201ll ). in this work 



the distance to maser galaxy NGC4258 was not used. 
Also, this work is based on a new, independent deter- 
mination of the distance to LMC, therefore, the cali- 
bration of Cepheid period-luminosity relation should 
be considered as independent one. For distance mea- 
surements differe nt d ata were also used (see details in 



Rics s et al.ll201ll and Frccdm an et al.l 12012). There 



fore, apparently, the Hubble constant measurements 
presented in these two works should also be consid- 
ered as independent ones. 

These two measurements combined in assumption of 
their independence and gaussian errors give the value 
Hq = 74.1 ± 1.8 km s _1 Mpc _1 . This measurement is 
designated below as Ho- We note, ho wever, that new 



Hubble constant measurement from Freedman et al 



(|2012l ) do not produce strong changes in the con- 
straints presented below. The constraints obtained 



1 http: / /hea. iki.rssi.ru/ 400d / cosm / 



using Hq measurement from IRiess et al.l (|201ll ) only 
are also given below for reference. 

Baryon acoustic oscillations 

The results of new, considerably improved measure- 
ments of baryon acoustic oscillations made using the 

2 http: / /lambda. gsfc.nasa.gov/ 
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data of large spectroscopic surveys of galaxies, pub- 
lished recently, are also used in our work. They 
inc lude the reprocessed data of SDSS Data Release 
7 (IPadmanabhan et all 2012), the data of Wi gqleZ 
(jBlake et all 1201 lh and 6dF (jBeutler et all l201lh sur- 
veys, as well as the measurements made using the data 
of S DSS Data Relea s e 9 (B OSS survey, CMASS sam- 



Relativistic energy density in early Universe is usu- 
ally parametrized using the effective number of neu- 
trino species, N e ff\ 



pie, lAnderson et all 12012). All these data taken to- 



gether are designated below as BAO. 

ACDM MODEL 

The constraints on the mean matter density in Uni- 
verse, f2 m , the linear density fluctuations amplitude, 
as, and Hubble constant, H$, which are obtained 
using various cosmological d atasets in standard six- 



parameter ACDM model fe.g-. lLarson et all 1201 lh are 
shown in Fig. [TJ One can see that there is a dis- 
crepancy between the direct Hubble constant measure- 
ments and the measurement obtained from BAO data 
in assumption of ACDM model (see left panel of the 
Figure). One can also see that new BAO data pro- 
duce a discrepancy between the measurements of as 
from CMB+BAO and from galaxy cluster mass func- 
tion data (see right panel of Fig. [1]). 

The discrepancy in Hubble constant measurements 
appears most prominently in BAO data from SPS S 
DR9 (see Fig. 21 and Fig. 30 in lAnderson et al.ll2012h . 
This d iscrepancy was noticed and was discussed in de- 
tail by lAnderson et aL ( 2012 ). In this work BAO data 
are c ompared to t he me asurement of Hubble constant 
from iRiess et al. ( 2011 ). When anot her, more recent 



measurement of Hubble constant from Freedman et al 



(120121 ). is taken in account, statistical significance of 
this discrepancy becomes more significant. It is this 
discrepancy, which is shown in the left panel of Fig. [TJ 

ADDITIONAL NEUTRINO SPECIES 

The underestimated measurement of Hubble con- 
stant can be obtained from BAO data in assumption 
of ACDM cosmological model with standard number 
of neutrino species, if in real Universe the extra com- 
ponent in relativistic energy density is present at early 
epoch in addition to photons and three known neutri- 
nos. In this case the Universe expansion rate is in- 
creased during the radiation-dominated era, and the 
size of sound horizon is therefore decreased. If one use 
ACDM model with standard number of neutrinos, the 
size of sound horizon would be overestimated, as com- 
pared to its real size. Since BAO observations give the 
measurement of the distance in units of sound horizon, 
Hubble constant, measured in this way will be under- 
estimated. Useful discussion on relat ed su bjects can 
be found in Eisenstein fc White! (|2004l ) and lHou et aT 
(1201 lh . 



Pr 



7 / 4 \ 4/3 

+ 8 n' 



where p 7 — energy density of photons. The con- 
straints on D, m , ag and Ho in ACDM model with free 
effective number of neutrino species added are shown 
in Fig. [2J From this Figure one can see that the dis- 
crepancies in Hubble constant measurements are elim- 
inated in this cosmological model, as expected. 

ADDITIONAL NEUTRINO SPECIES AND 
NON-ZERO TOTAL NEUTRINO MASS 

In Fig. [5] one can also see that even if additional 
neutrino species are introduced into ACDM model, not 
all the discrepancies are eliminated. From the right 
panel of this Figure one can see that the combined data 
on CMB power spectrum, BAO and Hubble constant 
measurements suggest significantly higher value of as, 
as compared to the cluster mass function data at the 
same value of f2 m . 

The measurement of as from CMB+BAO+Hq data 
is inferred from the amplitude of CMB temperature 
fluctuation, which contain the information on density 
fluctuation amplitude at high redshifts, z f» 1000. On 
the other hand, the measurement of as from the galaxy 
cluster mass function data reflects the value of this 
quantity at recent epoch. Therefore, this discrepancy 
in as measurements can be interpreted as a result of 
a suppression of density fluctuations due to non-zero 
total neutrino mass. 

The constraints on f2 TO , as and Hq in 
ACDM+Nes+Emv model arc shown in Fig. [3] 
One can see that in this case all cosmological datasets 
are completely consistent with each other. The 
consistency is achieved at the expense of the de- 
tection of larger then standard effective number of 
neutrino species and non-zero total neutrino mass. 
The constraints on total neutrino mass and effective 
number of neutrino species obtained with all the data 
combined are shown in Fig. [4] and [5] 

If all considered cosmological data are used 
(CMB+BAO+Hq+CL), the change of \ 2 when two 
parameters, Sm^ and N c s, are introduced into the 
cosmological model turns out to be A\ 2 = 13.0, which 
corresponds to ~ 3.2a significance. In this case the 
following measurements of effective number of neu- 
trino species and total neutrino mass are obtained: 
N e $ = 4.03 ± 0.36, Y,m v = 0.49 ± 0.17 cV. Here 
and everywhere below the systematic uncertainties 
of galaxy cluster mass function measurements from 
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ACDM ACDM 




Fig. 1. The constraints on Q m , and Ha in ACDM model, obtained using various cosmological datasets. The constraints from 
CMB+CL data are shown with solid lines, from CMB+Hq data — with dashed lines, from CMB+BAO data — with dotted lines. 
Shaded regions show model independent constraints (at lc level) from direct measurements of Hubble constant (left) and for the 
combination of trg and Q m from galaxy cluster mass function (right). 
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Fig. 4. The constraints on Em„ and N e g in ACDM+N e ff+'Em v model. The dashed contours show the constraints ob- 
tained using CMB+BAO+Hq data, solid contours — the constraints with galaxy cluster mass function data taken in account, 
CMB+BAO+Hq+CL. In the left panel the contours obtained with the same data assuming that the masses of clusters in 
IVikhlinin et al.l j2009al ) are underestimated by 9% are also shown with dotted lines. 
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Fig. 5. The constraints on Em„ and N e g in 

KCDM+N c g +T 1 m l , model. The dashed contours show 
the constraints obtained using CMB+BAO+Hq data, solid 
contours — the constraints with galaxy cluster mass function 
data taken in account, CMB+BAO+H +CL. 



Fig. 6. The constraints on Em„ and N e s in 

ACDM+N c ff +Em„ model in the case when only one 
neutrino type is massive. The contours are the same as in 
Fig. [5] 



Vikhlinin et al.l (|2009bf ) are included in the error of to- 
tal neutrino mass measurement. Note, that the mea- 
surement of effective number of neutrino species do 
not depend on cluster mass function data (see right 
panel in Fig. @|. For the ca se, when the mea sure- 
ment of Hubble constant from Riess et al. | (1201 lh only 
is used, the change of \ 2 when two parameters are 
added to cosmological model is A\ 2 = 11-9, which 
corresponds to ps 3. Oct significance. In this case the fol- 
lowing measurements are obtained: N e g = 3.89 ±0.39, 
Sm„ = 0.47 ±0.16 eV. 

When these constraints are obtained, the total neu- 
trino mass is assumed to be distributed in equal parts 
among three neutrino types (this is the default setting 
in used version of CosmoMC software). Different distri- 
bution of neutrino mass can produce notable changes 
in total neutrino mass constraints since, with the same 
total neutrino mass, more massive neutrinos become 
non-relativistic at earlier time. Correspondent differ- 
ences in total neutrino mass constraints are indeed no- 
table when current obser vational data are used (e.g., 
Burenin. Vikhlininl . 120121 ) . For this reason we also ob- 
tained constraints in the case when only one neutrino 
type is massive. These constraints are shown in Fig. El 
In this case the upper limit for total neutrino mass is 
somewhat relaxed, while the statistical significance of 
the measurement of non-zero neutrino mass remains 
to be approximately the same. 

In this case the change of \ 2 when two parameters 
arc added to the model is A% 2 = 14.6, which corre- 
sponds to w 3.4ct significance. In this case the follow- 



ing measurements are obtained: N c g = 4.32 ± 0.39, 
Era,; = 0.60 ± 0.22 eV. If the H ubble constant mea- 
surement from Riess et al. | (l201lh only is used, we get 
A% 2 = 12.8 which corresponds to sa 3.1ct significance, 
and in this case the measurements: N e g = 4.20 ±0.46, 
Em„ = 0.58 ± 0.23 eV are obtained. 

DISCUSSION 

As it is shown in Fig. [TJ all the discrepancies in 
cosmological data in ACDM model with zero neutrino 
mass and standard number of neutrino species appear 
mainly due to the usage of new baryon acoustics oscil- 
lations data. In fact just these data require the intro- 
duction of additional parameters into ACDM model, 
like non-zero total neutrino mass and larger than stan- 
dard number of neutrino species. If these BAO data 
are excluded from the consideration, all other data re- 
main to be consistent with standard ACDM model (see 
Fig. P. 

We note that in addition to the data used in our 
work, the other BAO observations exist where the 
lower value of Hubble constant is obtained in ACDM 
model, as compared to direct Hubble constant mea- 
surements. For example, similar discrepancy was 
found in the measurement of angular diameter dis- 
tance from BAO observed using photomet ric redshifts 
of luminous red galaxies (|Seo et all 120121) . The data 
on BAO observations in the transmitted flux fraction 
in the Lya forest of h igh redshift quasar s obtained 
using SDSS DR9 data (jBusca et all l2012h published 
recently, also give underestimated Hubble constant 
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in ACDM model (see. Fig. 21 in iBusca et all 120121) . 
Note, that these BAO data are obtained using signif- 
icantly different methods, as compared to BAO mea- 
surements in large spectroscopic galaxy surveys used 
in our work. 

As it was discussed above, this discrepancy was no- 
ticed earlier, and it was also found that it can be 
eliminated if extra energy density of relativisti c mat- 
ter is present in early Universe dMehta et al. . 2012 



Anderson et all 120121 : iFreedman et all 120121 ). Inde 
pendent indications for the possibility of larger than 
standard value of N e g, were also obtained earlier using 
the data on CMB anisotropy at smaller a ngular scales 
( Dunklev et al. 2011 ; Keisler et al]2011, see also dis 



cussion in Hamann et al. 20 111 Joudaki et al. 2012 : 
Riemer-Sorensen et al.ll2013T ). In this case the addi- 
tional radiation energy density in early Universe is ob- 
served using different physical effect, by the measure- 
ment of reduced power in CMB power spectrum damp- 
ing tail due to larger photon diffusion angular scale. 
These measurements are degenerate with primordial 
helium abund ance and running spec tral i ndex (see de 



tails, e.g., in iDunkley et al.l 120111 and Keisler et al 



20111 ). The data from lKeisler et al.l (|201ll ) are used in 
our work as well, and they also have some influence on 
our N g measurement, however, this influence is not 
dominant. As it was discussed above, the measure- 
ment of A^off, presented in our work, is based mainly 
on the new data on BAO observations compared to 
the direct Hq measurements. 

From Fig. Q] one can also see that current data on di- 
rect Hubble constant measurements and cluster mass 
function are consistent with each other in standard 
ACDM model. In order to bring these data in accor- 
dance to recent BAO observations one need to decrease 
the measured value of Hubble constant by 7-10% and 
to simultaneously increase the value of as measured at 
the same VL m also by 7-10%. 

In our work we use two independent measurement 
of Hubble c onstant which are in good agreement wit h 
each other ( Riess et all 1201 It IFreedman et al. 1. 120121) . 
Each of them is accurate to approximately 3%, in- 
cluding systematic errors. These are the data which 
give the measurement of larger than standard neutrino 
species when are compared to BAO observations. 

The data on galaxy cluster mass function in their 
turn, when compared to the data on BAO observa- 
tions, give the measurement of non-zero total neutrino 
mass. These constraints are affected by systematic er- 
rors of cluster mass function measurements. However, 
these data would be consistent with zero total neutrino 



mass o nly if the masses of clusters in IVikhlinin et al 



(2009a) were underestimated by approximately 30%. 



d ata (CMB+BAO+Hp+ CL), in case if cluster masses 
in Vikhlinin et al. ( 2009a|) are underestimated by 9%, 



the adopted systematic error for cluster mass scale cal- 
ibration. 

Systematic errors in determination of cluster mass 
scale is one of the main systematic uncertain- 
ties in cluster mass function measurements. In 



Vikhlinin et al.l (|2009bl ) the measurements of cluster 
masses were based on the temperature and mass of hot 
intracluster gas and were calibrated using hydrostatic 
mass measurements. Systematic error of cluster mass 
scale calibration was estimated as 6M/M « 0.09, us- 
ing the comparison of hydrostatic cluster masses with 
those bas e d on w eak lensing measuremen ts, taken from 



Hoekstral (|2007l) and lZhang et al.l (|2008l) . More recent 



works, where the additional weak lensing data are pre- 
sented, are in general agreement with approximately 
10% systematic error for cluster m asses measured us- 
ing the data of X-ray observations (llsrael et al. 
20121: iMahdavi et all l2012t lAppleeate et all 120121) . 



CONCLUSIONS 

In our work we show that the discrepancies between 
different cosmological datasets in the determination of 
Hubble constant and in measurements of density fluc- 
tuation amplitude in assumption of standard ACDM 
cosmological model, can be eliminated if the additional 
neutrino species and non-zero total neutrino mass are 
introduced into the cosmological model. In this case, 
the discrepancy in the determination of distance scale 
between baryon acoustic oscillations data and direct 
Hubble constant measurement is eliminated by the 
assumption of larger than standard effective number 
of neutrin o species, which was already noticed earlier 



Csee. e.g.. iMehta et al l. l2012t lAnderson et all l2012t 
Freedman et all l2012t ISeo et all l2012t) . The remain- 



One can see this from the left panel in Fig. @J where 
the dotted lines show the contours for all considered 



ing discrepancy in as measurements between combined 
CMB, BAO and Hubble constant data and the data 
on galaxy cluster mass function is eliminated by the 
assumption of non-zero neutrino mass. 

The change of x 2 when two parameters, N c g and 
£m„, are introduced into the cosmological model cor- 
responds to w 3cr significance level. The model with 
approximately one additional neutrino type, N e g ?s 4, 
and with non-zero total neutrino mass £m„ 0.5 pro- 
vide the best fit to the data. In model with only one 
massive type of neutrino the upper limits on neutrino 
mass are slightly relaxed. 

We emphasize, that ACDM model with standard 
number of neutrino species and zero neutrino mass 
appears to be no longer consistent with cosmological 
data due to the usage of the data on baryon acous- 
tic oscillations, published recently. In future these 
results may be independently confirmed using signif- 
icantly improved data on CMB anisotropy at small 
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Fig. 7. The constraints on £m„ and N e g in 

KCDM+N eS +T,m u model, obtained using WMAP9 data 
and also with the data on CMB lensing potential included. 
The dashed contours show the constraints obtained using 
WMAP9+SPT+CMBLcns+BAO+H data, solid contours - 
the constraints with galaxy cluster mass function data taken in 
account, WMAP9+SPT+CMBLens+BAO+H + CL. 



Fig. 8. The constraints on Em„ and N e g in 

ACDM+N c ff +Tim u model using all considered data, ex- 
cluding the data from SPT. The dashed contours show the 
constraints obtained using WMAPT+BAO+Hq data, solid 
contours — the constraints with galaxy cluster mass function 
data taken in account, WMAP7+BAO+Hq+CL. 



angular scales. In very near future the measurements 
made with Planck space observatory should be pub- 
lished which will probably clarify this issue. 

It is known that the existence of additional light 
neutrino species, which could produce oscillations 
with Am 2 ~ 1 eV 2 , was suggested to explain the 
results of the experiments of neutri no oscillations 
searches at short baselines — LSDN (Aguilar et al 



200lh . MiniBooNE (jAguilar-Arevalo et all l2010h 



Mention et all 201 lh. a 


(lAbdurashitov et a 


.. 20061: 


Giunti & Lavedeil 


20111). 1 



re- 



2011 



gallium a n omal; 



201 



In order to explain the 
results of these experiments the possibility of the 
existence of one or two lig ht sterile neutrino wa s 
discussed (see details in, e.g.. lAbazaiian et all 120121) . 

These neutrinos can be constrained from the cos- 
mological data as well. The constraints presented in 
our work are consistent with the existence of one ad- 
ditional type of neutrino. Since in this case Am ~ m, 
the masses of all other neutrinos should be significantly 
lower, and to constrain these neutrinos one should use 
the model with only one massive type of neutrino. 
These constraints were obtained above and they are 
as follows: N cfi = 4.32 ± 0.39, £m„ = 0.60 ± 0.22 eV. 

JVote 

When this paper was already submitted for publi- 
cation, few papers were published, where some new 



constraints on number of neutrino species and to- 
tal neutrino mass based on new cosmological data 
were presented. They are the papers on cosmo- 
logical constraintsobtained from nine year WMAP 
data (Hinshaw et al ], 120121) and from new measure- 
ments of CMB anisotropy at high multi poles using 
the d ata from South Pole Telescope ( SPT AKou et all 



2012 l) and Atacam a Cosmology Telescope (ACT, 



Sievers et all 12013 ) 



In their paper on nine year 



WMAP data, 
Hinshaw et all (|2012h . obtained N eS = 3.84 ± 0.40 in 
model with free N e s and zero total neutrino mass, 
in agreement with our constraints. Note, that in 
this work CMB data include also the data on power 
spectrum of CM B gr avitational lensing potent ial from 



Das et al.l (|2011l) and Ivan Engelen et al . (2012). which 



contain the information on the density fluctuation am- 
plitude at redshifts approximately 0.5 < z < 5. If 
the data on density fluctuations amplitude at recent 
epoch are added to other current cosmological data 
in ACDM+A r c fj model with zero neutrino mass, the 
lower value of N e g is measured due to the degener- 
acy between N e g and Em ^, which appears in this case 
(Bur enin. Vikhlininl [20121 see also, e.g., contours for 
Neff-Hm^ in Figures above). 

The effect of these new data on the constraints 
presented above is shown in Fig. [7J In this Fig- 
ure the constraints on number of neutrino species 
and total neutrino mass obtai ned using the la t est 9 - 
year WMAP data (WMAP9, IHinshaw et all l2012h 
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and also the data on CMB l ensing pow e r spec trum 



amplitude ( CMBLens) , from iDas et al.l (|2011l ) and 



van Engelen et al.l (|2012T ) 



which was taken in account 
as a gaussia n prior, C2lm = (3 . 17 ± 0.45) x 10 -8 , as it 
was done bv ICalabrese et all ( 20131 ). are shown. The 
new data only slightly change the results presented 
above. In this case the significance of the inclusion 
of two parameters, N e s and T,m v in the model is 
~ 3.3a, and we obtained the following measurements: 
N cS = 4.02 ± 0.35 and Sm„ = 0.50 ± 0.15 cV. If 
only one type of neutrino is massive, the correspon- 
dent significance is approximately the same and we 
obtain: N cS = 4.12 ± 0.36 and Ym v = 0.58 ± 0.21 eV. 
Note, that when CMB lensing power spectrum data 
are used, non-zero neutrino mass is obtained with 
lower significance even without clusters mass function 
data (dashed contours in Fig. [7]). 

The cons traints based on n ew data from South Pole 
Telescope (|Hou et all . l2012h in ACDM+7V cff+ £ mi/ 
model are generally consistent wit h the resul t s, ob - 
tained in our work. From Fig. 18 of iHou et al.1 (|2012j ) 
one can see that, similarly to our work, the mea- 
surement of larger than standard effective number of 
neutrino species is based mainly on the new data on 
BAO observatio ns (see also the discussion in §10 of 
Hou et al.l 120121 ) . Also, as in our work, the measure- 



ment of non-zero neutrino mass comes from the data 
on galaxy clusters m ass functio n . The refore, as com- 
pared to our work, IHou et al.l (|2012l ) used approxi- 



mately similar cosmological data and obtained approx- 
imately similar constraints on neutrino properties, as 
expected. 

In their work, based on the new data of At- 
acama Cosmology Telescope (ACT), ISievers et al 



( 2013T ) did not find additional neutrino types when 
both ACDM+N c s and A CDM+N e s+Em u models 
are considered. Note, that not complete BAO dataset 
is use d in this work (o nly th e data from SPSS 
DR7, IPercival et ail l201Ct 6dF, iBeutler et all 12011 



and SPSS DR 9, Anderson et al.l 20121 ). Also, as in 

11- 



Hinshaw et"aH <l2012h and IHou et all hoitfl . the H 
measurement from lFreedman et al.l ( 20121 ) is not used. 
Apparently, in the new data from ACT there are 
some statistically insignificant discrepancies with SPT 
data which bring the measurement of N e s to some 
lower value. However, as it was discussed above, in 
our work the measurements of larger than standard 
number of neutrinos and non-zero neutrino mass are 
based mainly on recent BAO data and do not depend 
strongly on the data on high-/ CMB measurements. 

In order to better show this, in Fig. [8] the con- 
straints on N e g and Em„ are given, which were ob- 
tained using all the data considered in our work, ex- 
cluding the high-Z CMB data from SPT (i.e., using 
the WMAP7+BAO+H +CL dataset). There is no 



much change in the significance of the detection of 
N e g > 3 and Y>m v > 0, the main effect of the ex- 
clusion of SPT data is the relaxed upper limit for 
N e ff. In this case, the change of x 2 when two param- 
eters, N c g and Em„, are introduced into the model 
is A\' 2 = 11.1, which corresponds to « 2.9er signifi- 
cance, and the following measurements are obtained: 
N oB = 4.44 ± 0.50, Sm„ = 0.49 ± 0.18 eV. In the 
case, when only one neutrino type is massive, A% 2 
have the same value, and the following measurements 
are obtained:^ = 4.55±0.50, Y>m v = 0.54±0.24 eV. 
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